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Reply: In our Letter [l| we presented a numerically 
exact, polynomial-time quantum computer algorithm for 
simulating the class of general pairing Hamiltonians (with 
arbitrary coupling constants U^ ; ) 
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[see Eq. (1) in l] for definitions]. The Comment by 
Dukelsky et al. |2| reviews an interesting but numerically 
approximate classical algorithm [density matrix renor- 
malization group (DMRG)] which they claim to be appli- 
cable for the same problem. It is this difference between 
numerically approximate and exact algorithms that es- 
sentially renders irrelevant the criticism of our Letter 
in the Comment Q. There is no known classical algo- 
rithm for simulating the class of general pairing models 
that is both efficient (polynomial) and numerically exact. 
There is such a quantum algorithm: the one proposed in 
our Letter. 

Unfortunately there is a misleading and confusing 
statement in the Comment 2]: Dukelsky et al. write 
in the third paragraph that (all emphasis ours) "It was 
shown in this context that only by resorting to the exact 
numerical solution [4], using the DMRG However, 
the DMRG is not an exact algorithm. It is an approx- 
imation, as in fact stated by the authors in their Com- 
ment a few lines after their claim that the DMRG is ex- 
act. Indeed, citing Dukelsky and Sierra's own paper |3| 
(their Ref. [4], Ref. [14] in our Letter): "We show that 
the DMRG gives an accurate approximation to the exact 
ground state if the block is taken to be the set of particles 
while the environment is taken to be the set of holes." 
And on p. 174 of the same paper they write: "The DMRG 
is a variational method and in the region under study we 
expect our results to coincide with the exact ones with a 
relative error less than 10 -4 ." 

These statements by the authors of the Comment 
clearly demonstrate that the DMRG method is not a 
numerically exact solution of the pairing Hamiltonians. 
Hence the term "exact" in the the above quote from their 
Comment is not correct, and the basis for their criticism 
of our Letter is invalidated. 

The authors have provided evidence that the DMRG 
can give a good approximation to the low-lying energy 
spectrum of the reduced BCS Hamiltonian (with constant 
pairing interaction = V r ), and can do so in O(N) 
steps We do not dispute this, and made no state- 
ment to the contrary in our Letter, since we were only 
concerned with a numerically exact algorithm. However, 
in spite of their claim to the contrary in their Comment 
("...the DMRG can easily accomadate arbitrary pairing 
matrix elements."), there is presently no evidence to sup- 
port the hope that the quality of the approximation and 
the efficiency of the algorithm will remain as good for 
general pairing Hamiltonians (i.e., VL^ ^ V r for all m,l). 


The burden of proof is upon the authors of the Comment. 
Moreover, even if they could demonstrate this, the basic 
point remains that theirs is a numerically approximate al- 
gorithm while ours is exact. Indeed, as we stated in our 
Letter, for half-filled states the dimension of the Hilbcrt 
space grows as N\/[(N/2)\] 2 (for even N), which is super- 
exponential in N, and there is no known way to avoid 
this divergence in a classical algorithm that is numeri- 
cally exact. In contrast, we have given a quantum com- 
puter algorithm that is numerically exact (in the sense 
that it exactly diagonalizes the pairing Hamiltonian), re- 
quiring only 0(N 4 ) steps, for arbitrary N. Any approx- 
imate method, including DMRG and projected BCS U 
must eventually be checked against numerically exact al- 
gorithms if and when we have quantum computers. 

In the second part of the Comment, the authors con- 
sider a pairing model that includes a monopole interac- 
tion [third term in their Eq. (1)], not included in our 
pairing Hamiltonian, so it does not apply to our Let- 
ter. They point out that if the coupling constants Vl 
and Vjj (in their notation) satisfy certain relations then 
the model is exactly solvable. This is an elegant result 
of which we were well aware, but that bears no relation 
to our work, which deals, of course, with the cases for 
which no analytical solution is possible. However, Dukel- 
sky et al. write 0: "We believe that most of the physi- 
cal problems can be modeled with a pairing Hamiltonian 
within the integrable subset ..." , indicating that the cases 
that are not analytically solvable are somehow typically 
not physicallyrelevant. This conclusion is invalidated by 
their own parameter count: the ratio of integrable to 
non-integrable models is [(6N + 3)/(2N 2 - N)} — >• as 
N — > oo, so the class of integrable models occupies a 
negligibly small fraction of the parameter space in the 
interesting regime of N ^> 1. 

In conclusion, the Comment discusses approximate 
classical algorithms and analytically solvable instances 
of the pairing model, which are certainly important, 
but are not relevant in the comparison to a general and 
numerically exact algorithm for this problem, such as 
ours. 
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